Initial Decoherence and Disentanglement of Open Two-Qubit Systems 
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Abstract 

In this paper we investigate a open two-qubit model whose dynamics is not exactly solvable. When the initial state is 
the maximum entangled state, as the exactly solvable open two-qubit model [D. Tolkunov and V. Privman, Phys. Rev. A 
71, 060308(R) (2005)], the decay of entanglement of formation of the model, expressed by concurrence is also governed by 
the product of suppression factors describing decoherence of the subsystems (qubits). However, if the initial state is not the 
maximum entangled state, its concurrence will decrease faster than the product of the suppression factors describing decoherence 
of the qubits. 
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I. INTRODUCTION 



Coherence and entanglement are two remarkable fea- 
tures of quantum systems, and they are also dominating 
reasons why quantum computation and quantum com- 
munication differ from the classical ones 1]. For exam- 
ple, quantum coherence can lead to natural parallel com- 
putations which can enhance efficiency for solving some 
complex problems by using effective quantum algorithms. 
Quantum bit (qubit) is a key block for building quantum 
computers. However, the interactions of qubits with their 
environment will induce loss of the coherence, decoher- 
ence. Decoherence is considered a central impediment to 
fabricate quantum computers. In quantum computation 
and quantum communication, the entanglement, nonlo- 
cal correlation between the quantum subsystems is also 
needed. The entanglement is recognized as an important 
physical resource in quantum information transmission 
and processing. Many protocols in quantum communi- 
cation and quantum computation are based on the en- 
tangled states. It is also shown that couplings of the 
quantum systems and its subsystems to their environ- 
ment will result in the loss of the entanglement, and this 
loss cannot be restored by local operations and classi- 
cal communications 0, ll| • Both the decoherence and the 
loss of the entanglement may result from the interactions 
of the quantum systems with their environment. Then, 
what relations are there between the two quantities? In 
this paper we will investigate the relationship of the de- 
coherence and the loss of entanglement for a non-exactly 
solvable two-qubit model. As Ref. 07 Hi , let us refer to 
two subsystems, S**^^^ and S**^^-*, of the combined system, 
S. The evolutions of the coherence of the subsystems 5*^^-' 
and 5^^^ in their bathes can be described by suppression 
factors, < ^'^'^-^ [t) < 1. On the other hand, the evo- 
lutions of the entanglement between the subsystems S*^^-* 
and 5'(2) can be described by entanglement of formation 
expressed by concurrence 'l5"|. In large times the de- 
cay of the coherent terms are exponential, then the decay 
rate of the whole system is the summation of the rates of 
the subsystems 0,11 [l||. Yu and Eberly ij found that 



for a exactly solvable model in which the decoherence is 
caused by pure dephasing process, the concurrence de- 
cays faster than the quantum dephasing of an individual 
qubit. A more physical model, two entangled atoms in 
pure vacuum noise, is investigated recently by the same 
authors They found that the disentanglement time 
is shorter than the usual spontaneous lifetime. V. Priv- 
man et al. 0, 13 investigated another dynamically solv- 
able model. They shown that the decay of concurrence of 
the model is governed by the product of suppression fac- 
tors describing decoherence of the subsystems (qubits) 
in a short time. It is interesting that whether or not 
the relationship of the concurrence and the suppression 
factors describing the decoherence can be held for other 
open quantum systems. In this paper we shall investigate 
a different open two-qubit model whose dynamics is not 
exactly solvable. By using a short-time approximation 
we can obtain the evolution of reduced density matrix of 
the system. It will be shown that as the initial state of 
the two-qubit system is the maximum entangled states 
the decay of concurrence, namely the decay of the de- 
gree of entanglement of formation can also be governed 
by the product of suppression factors describing deco- 
herence of the subsystems (qubits) in a short time. If 
the initial state is not the maximum entangled state, the 
concurrence will decrease faster than the product of the 
suppression factors describing the coherence of the qubits 
does. 



II. DECOHERENCE AND THE LOSS OF EN- 
TANGLEMENT 



Suppose the open two-qubit system has Hamiltonian 



(1) 
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where 

k k 

ffj = alY,(gl*hl+glhl^). (2) 

fe 

Here, we use the subscripts or the superscripts r — 1,2 
to label the qubits. iJJ and are the Hamiltonian of 
qubits and their bosonic bathes and i/J are the inter- 
actions of the qubits with their bathes 11]. Where we 
assume that each qubit interacts with its own bath. This 
is not a exactly solvable model. If we do not consider the 
interactions between the qubits the evolution operator of 
the combined system can be expressed as 

U = Ui® U2. (3) 

The evolution operator of the single qubit is 

where t — t/H. Due to non-conservation of Hg in this 
system, the evolution operator cannot be in a general 
way expressed as e~'^=*e~*^^^"'"''^s-''. But in the sort- 
time approximation, the operator can be approximately 
expressed as 0, 

Ur - e-''^^*/2e-'(^?+««)*e-*^»^*/2 + o{t^). (5) 

It has been proved that the expression is accurate enough 
as the time being short to the characteristic time 0| . So 
the elements of the reduced density matrix p (t) in the 
basis of operator Hs can be expressed as 

p;;,„ = Tr5(v.„Je-^='/2e-(«?+^«)*e-''^r*/2^(0) 

Here, we suppose the initial state of the system be 
(0) = p'' (0)(g)e'' where p'' (0) and 6'' are initial states 
of the qubit and its environment. The latter is the prod- 
uct of the bath modes density matrices 9\ . In the initial 
states, each bath mode k is assumed to be thermalized, 
namely, 

" Tr, (e-Z^A^J) ' 

where P — 1/kT, T is the temperature and k is the Boltz- 
mann constant. 



A. Decoherence and decay of coherence 

In this subsection we shall investigate the decoherence 
of each qubit due to the interaction of the qubit with 
its own environment, a bath. Here, we denote the sup- 
pression factors describing the decoherence of each qubits 



with . The relationship of the suppression factors 6^ 
with the usual term decoherence is = L {1 — S^) , 
where L is the initial coherent terms of the density ma- 
trix, namely, the off-diagonal elements of the density ma- 
trix of the initial state for the r ~ th qubit. L = 1/2 for 
the initial states with maximum coherent terms. We call 
L5^ the decay of coherence of the r — th qubit. Usually, 
the environment is assumed to be a large macroscopic 
the interaction with it leads to the thermal equilibrium 
at temperature T. In this case, Markovian type approxi- 
mations can be used to quantified the decoherent process 
and it usually yields the exponential decay of the density 
matrix elements in the energy basis of the Hamiltonian 
i/J. Thus, the decay rates will be additive. In this time 
scale the measures of entropy and the first entropy can 
be used for quantifying the decoherence. But the deco- 
herence of the qubit gate operations cannot be character- 
ized by this methods because the time of the elementary 
quantum gate operation is much shorter than the thermal 
relaxation time. It has been shown that the norms ||cr''||j;^ 
is useful for describing the decoherence of the short-time 
evolution Here is the deviation operator defined 
as 

it) = (t) - p[ it) , (8) 

where p'' {t) and (t) are density matrixes of the "real" 
evolution (with interaction) and the "ideal" one (without 
interaction) of the r — th qubit. We can use the norm 
||f7''|U describing the decoherence and the norm is defined 
as JJI 

For a qubit, the norm can be given by 

hlx^^Wlof + Wlif- (10) 

It is shown that for a given system, the norm ||cr''|l_^ in- 
crease with time, reflecting the decoherence of the sys- 
tem. However, in general it is oscillated at the system's 
internal frequency. Thus, the decohering effect of the 
bath is better quantified by the maximal operator norm 

D^{t)^ sup {\\a^{t,p^m\\,). (11) 

p'-(O) 

For our investigating model, by using Eqs.© and lllfl we 
can obtain the decoherence of the qubit as 

D^it) = l{l-e-^^^W), (12) 

where 

(i) ^ 2 E ^ sin^ ^ coth M . (13) 

It is shown that the suppression factor describing the 
decoherence of the open qubit r is 
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B. Loss of the entanglement and decay of the en- 
tanglement 

Numerous measures of entanglement have been con- 
sidered over the years. For quantum information con- 
tent, the entanglement of formation has been a widely 
accepted measure of entanglement. The measure can 
measure the degree of entanglement not only for pure 
states but also for mixed states. The concurrence is a 
quantity monotonically related to the entanglement of 
formation. For a pure or mixed state, p^, of two qubits, 
one can define the spin-flipped state. 



~Ps = i.<^y®'^y)P*s {<^y®<^y)^ 
and the Hermitian matrix, 



'PsP. 



srs\ rsi 



(15) 



(16) 



with eigenvalues Ai=i,2,3,4. Here, p* denotes the complex 
conjugation of p in the standard basis and Oy is one of 
the Pauli matrixes. The concurrence, C, is defined by 

m 



C — max < 0, 2 max 




From Eq.® and helped with operator-algebra techniques 
we can obtain the evolutions of the reduced density ma- 
trix as 



where 

Pll it) 

Poi W 



The concurrence of the initial state is C(0) = 
2 \a\ / + ■ At first, we investigate an especial 
case, namely, the case of a = 1, which denotes the ini- 
tial state is a maximum entangled state. We can easily 
obtain that after time t the density matrix of the open 
two-qubit system becomes 



Ps it) = 



1 



where A 



A Ae- 

B B 

B B 

Ae**^' A 

,-4Gi(t)-4G^(t) 



-itEj 



(23) 



and B 



e 4*5^ (*) 4G^(t)^ Then, we can obtain the eigenvalues of 
the product p^ (t) (t) as 



P 
and 



0. 



(24) 



(25) 



Finally, we can obtain the concurrence as 



(17) Cit) 



-4G^(t)-4G^(i) 



^pSo(l"e-^«^W)+ipii(l + e-^«'W). (19) 



C{0)5^ {t) (t) 
(26) 

We call the C (t) the decay of entanglement. The loss of 
the entanglement is 1 ~ C{t). If we set a = — 1 or ±i 
we can obtain some similar results. It is very interest- 
ing that in this case the suppression factors describing 
the decoherence of the qubits and the entanglement of 
formation between the two qubits still has the compact 
relationship as Ref.j^, j^]. But, if the relationship is still 
preserved for other initial states? In the following we 
shall investigate the cases that the initial states are not 
the maximum entangled states. When a ^ ±1 or ±z the 
operator calculations similar to above derivation becomes 
very complex, then we only be able to numerically inves- 
tigate this problem with a concrete bath model. Suppose 
the two bathes with Ohmic noise spectrum are the same 
for convenience, which do not loss the universality. The 
spectral densities of the bathes may be expressed as 



In the right hand side of Ea. H19|l we denote p[j (0) with 
p^o . If we do not consider the interactions between the 
qubits, the reduced density matrix of the combined sys- 
tem S becomes 



J (uj) = rjLoe 



J / u. 



It is well-known that when the summation in Ea. (|13|) is 
converted to the integration in the limit of infinite num- 
ber of the bath modes, one has 



p, (t) = pi (t) (g> p2 (t) 



(20) 



In the following, we set the system in a pure entangled 
state 



C (t) = 27] 

for the real g (uj) . 



2 ^t I3uj 
1 — coth — 



Ps (0) = |*) (*| 



at i = 0, where 



I*) = 



14^ lal 



[|01)+a|10)] 



(21) 



(22) 



Fig.l 
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In this paper we choose the cutoff frequency of the bath 
modes as oJc = 10^^ Hz. It is pointed that the di- 
mensionless strength of the dissipation is very small for 
some qubit environment. For example, for the Josephson 
charge qubit in Ohmic bath the value of 77 is about 10^^ 
|1(tiI17I |. In our calculations, if we take 77 = 10^6, the dif- 
ference of the decay of concurrence C [t) and the product 
of the suppression factors 5^ {t) 5^ (t) with a initial con- 
currence C(0) is very small. For convenience, where we 
denote S{t) := C(0)(5^ (t) 5^ (t) . When we increase 77 to 
10~^ we can clearly see that the decay of concurrence 
decreases faster than S{t) does. In Fig.l we plot C (t) 
versus time (with points), which is compared with S{t) 
(with lines) in different initial states where (a) a = 1, (b) 
a = 2, and (c) a = 3. It is shown that the concurrence 
is not equals to the product of suppression factors and 
its initial concurrence, namely S(t) except for a = I. 
In plotting of the figure we choose a time unit, ks for 
convenience. Where we set 1 ks = 1519.29 ps. From 
Fig.l we see that the durative time of our calculations 
is r = 12.15j)S, which is smaller than the characteristic 
time of the qubits So our calculations is accurate 

enough. 

It is shown that system 

2 

r=l 

where 

k 

Hi = ^:Y.{9rbl+glbl^), (27) 

k 



is also a non-exactly solvable model and it has the same 
dynamics to system Eas. (|ll2|l in short-time approxima- 
tion So the concurrence and the suppression fac- 
tors describing the decoherence of the system may have 
a similar relationship to system Eas. (|ll2|l . 



III. CONCLUSIONS 

In this paper we investigated a non-exactly solvable 
open two-qubit model. We obtained that in a short time 
the decay of entanglement of formation is governed by 
the product of the suppression factors describing deco- 
herence of the subsystems as the initial state of the two- 
qubit system is the maximum entangled state. This novel 
relationship is similar to the discover by V. Privman et 
at in 1^ i5| , where the open two-qubit system is exactly 
solvable. Our work shows that when the initial state is 
not a maximum entangled state, after the open two-qubit 
evolve a short time t the entanglement of formation, con- 
currence decrease faster than the product of the suppres- 
sion factors describing the decoherence of the two-qubit 
systern, which is similar to the discover by Yu and Eberly 
in |9|, [lOj . It is shown that when the dissipation is not 
very weak (the dimensionless strength of the dissipation 
1] ^ 10^^) the entanglement is distinctly more frangible 
than the coherence for the quantum systems. 
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IV. CAPTION OF FIG.l 

Fig.l: The concurrence C{a,t) (points) and S{a,t) = 
2a/(l + |a|^)5^ {t) 5^ {t) (lines) of a two-qubit system in 



their Ohmic bathes for different initial states (a) a = 1, 
(b) a = 2, (c) a = 3. Here we take dimensionless strength 
of the dissipation rj = 1 x 10~^, the cutoff frequencies of 
the bath modes ojc = 10^^ Hz. 
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